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Abstract 

We briefly review the 1 + 1-dimensional SU(N) gauge theory minimally coupled to an ad- 
joint scalar, based upon a dimensional reduction of 2 + 1-dimensional pure gauge theory, which 
approximates the dynamics of the transversely polarized gluon. The lightest glueball states are 
investigated non-perturbatively using light-front quantisation in the large-iV limit a la 't Hooft, 
and the physical meaning of the results elucidated. Comparison is made with recent lattice Monte 
Carlo results for 3-dimensional quenched QCD. 



1 Based upon talk given at the Spring Workshop on String Theory, Gauge Theory, and Quantum Gravity, Trieste, 5-7 
April 1995 
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1 Introduction 



It has been suggested by a number of authors |2], |3| that adjoint scalar fields coupled to 1 + 1- 
dimensional QCD provide a way of approximating in tractable form the features of boundstates 
arising from transversely polarized gluons in higher dimensional gauge theories^. Here we will 
re-examine in more detail using light-front quantisation the low-lying glueball states for the model 
describing the hr = gluons in pure 2+1-dimensional gauge theory^] in the large- iV limit, making 
comparison with recent preliminary results from lattice Monte Carlo simulation of 3-dimensional 
quenched QCD ||. Our aim is to assess the limitations of this approach in reproducing actual 
features of the higher dimensional theory. In the reduced model, apart from an expansion of 
the scale for ratios of boundstate masses, presumably due to level mixing as a result of broken 
rotational invariance, we find that the level ordering in the charge conjugation C = — 1 sector is a 
poor match with the lattice results, while the C = +1 sector fares much better and is qualitatively 
in agreement. We attribute the extreme precision of the valence gluon approximation in the model 
to transverse (physical) gluon interactions having to proceed via an instantaneous longitudinal 
intermediate gluon. 

Starting from SU(N) Yang-Mills theory in 2 + 1-dimensions 

S = ~ T2 I <? x Tri^F^ (1) 



4 <?3 ■> 

in the light-front gauge = one has in light-front co-ordinates x ± = (x° ± x l )/y/2 

S = ^ J dx 2 dx + dx-Tr (d+A 2 d-A 2 + ^(d-A + ) 2 + A+J + ^j (2) 

J± = i[A 2 , d-A 2 ]ij - d-d 2 A% ;i,j = l,...,N (3) 

The field A + is a constrained variable which does not propogate in light-front time x + , leaving 
only the transverse gluon A 2 as a physical degree of freedom. Much of the complexity of 3D gauge 
theory is due to the linear term in the longitudinal momentum current (||), so one might as a first 
approximation study the theory restricted to zero modes = only. This is equivalent to the 
1 + 1-dimensional adjoint gauge theory 



S R = J dx dx lr Tr 



(4) 



where g 2 = g 2 / J dx 2 , D a = d a +i[A a , .], and <\> = A 2 /g. A bare mass term has been added to (||) by 
hand in order to subtract a logarithmic ultraviolet divergence later in the calculation. Sr inherits 



2 Related work can be found in refs.Q]. 

3 The extension to 3+ 1-dimensional gauge theory with quarks will be the subject of a more extensive forthcoming 
publication. 
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a subgroup of the Poincare and gauge symmetries of the 2 + 1-dimensional theory. The residual 
internal symmetries of Sr resulting from gauge symmetries of 5 are just the 1 + 1-dimensional 
gauge symmetries 

A a -» UA a U ] + \{d a U)U ] , <f> -> UcpU ] (5) 

5^ acquires the 50(1, 1) subgroup of 50(2, 1) Lorentz symmetries of 5 and also has Z 2 symmetries: 
charge conjugation C (Afj — ► —Aj { ) induces a discrete symmetry (fiij — > — 0jj on the 1 + 1-scalar; 
there also remains the one-dimensional parity Pi symmetry (x 1 — > — x 1 ). We will in addition 
make use of other dynamical information in trying to match the 1 + 1-dimensional eigentates 
forming represenations of the above discrete symmetries with the representations J of the full 
spatial rotation group 50(2) in 2 + 1 dimensions. 

2 Boundstate Equations. 

When analysing the boundstates of Sr by light-front quantisation, taking the large-iV limit re- 
sults in additional simplifications, the calculation being similar to, but much richer than, the one 
performed by 't Hooft for quarks in the fundamental representation Q. The boundstate problem 
can be expressed entirely in terms of the Fourier modes aij(k + ) of the transverse gluon 4>ij{x~) 
(further details are given in refs.[Q, |3|). The physical Hilbert space is the light-front Fock space 
constructed from the creation operators a^(/c + ). Due to confinement and the large- N limit the 
light-front wavefunctions one needs to consider are the subset of singlet states of the form 

MP + ) >= £ f + dk+... dk+S (±k+- pA fn(kt N - n/2 ' K ] ^[a\kt) • • • aHtiW > (6) 

n=2 \m=l / 

Such states are eigenstates of P + , so the solution of the boundstate problem is reduced to the 
diagonalisation of P~ in this basis, the dispersion relation being M 2 ^ = 2P + P~ty . The coefficients 
f n of n-gluon states, cyclically symmetric in their arguments, satisy an infinite set of coupled 
Bethe-Salpeter integral equations 0, |2l M 

g z N g z N x\ A^fx\X2 

rx\+X2 

+ / dy {E[x 1 ,x 2 ,y]fn(x 1 ,x 2 , ,...,x n ) 





+ (D[x 1 ,x 2 ,y] - E[x 1 ,x 2 ,y])f n (y,x 1 + x 2 - y,x 3 , . . .,%)} 

rxi rxi—y 

+ / dy dz {F[xi,y,z]f n+2 (y,z,x 1 - y - z,x 2 , . . . ,x n ) 

Jo Jo 

+F[x 3 , -x 2 , -xi]/ n _ 2 (xi + x 2 + x 3 , x 4 , . . . , x n )} 

+cyclic permutations of (x\, x 2 , . . . ,x n ) (7) 
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D[xi x 2 y] = — ~ Xl ^ Xl +X2 ~ 2y ^ (8) 

4(xi + x 2 ) 2 y/xix 2 y(x 1 + x 2 + y) 

E r Xi x<2 y\ _ (xi +y)(x 1 + 2x 2 - y) ^ 

4(xi - y) 2 y/xix 2 y(x 1 + x 2 - y) 

F[x u y,z] = 1 gx l + y){x l -y-2z) | (gg - y - z){y - z) \ 

4Vxiyz(xi - y- z) \ (xi - y) 2 (y + z) 2 / 

where the longitudinal momentum fractions x m = k^ n /P + have been introduced. To these equa- 
tions we must add a number of provisos. As with the 't Hooft model || there are some bare mass 
ambiguities, since <j) has logarithmic and linearly divergent self-energies in 1 + 1 dimensions. The 
mass of an isolated gluon in higher dimensions does not get perturbatively renormalised as a con- 
sequence of gauge invariance, but is presumably pushed to infinity by non-perturbative confining 
effects. In 1 + 1 dimensions on the other hand colour flux is forced into tubes by the restricted 
space and confinement occurs perturbatively. The transverse gluon field 4> has divergent pertur- 
bative mass renormalisations. In singlet combinations of these gluons however, linear divergences 
cancel, resulting in the principal value nature of the (Coulomb potential) E integral in (^). For 
such 1 + 1-scalar particles the logarithmic divergences do not cancel and one must decide what to 
do with them. Adding a bare mass mo to Sr, which does not violate the reduced gauge invari- 
ance, we will adopt the prescription of ref.|| setting the renormalised mass m = in (Q). This 
is clearly not the only prescription one could take, but is perhaps intuitively closest to the three 
dimensional theory, which one is trying to model as accurately as 1 + 1-dimensions allow, since 
the spectrum is unbounded for m < while for m — > + a positive spectrum is given in units of 
the only dimensionful scale g 2 . More generally one might keep it as a free parameter and fit to 
a known spectrum in the spirit of ref.fij. Our choice of a perturbative (light-front) vacuum |0 > 
in any case implies that we are neglecting zero modes in this two-dimensional model — we see 
no reason or evidence for zero modes condensing in this case — the perturbative vacuum being 
sufficient for confinement in 1 + 1 dimensions. 

Simple approximate analytic solutions may be found to the equations (|7|) while numerical 
solutions may be easily obtained to any desired accuracy with the help of Mathematica and 
a workstation. By cutting the interval of allowed momentum fractions < x < 1 into x E 
(1/K,3/K,5/K, ■ ■ ■) for some integer K, the problem becomes one of finite matrix diagonalisation 
Q, the continuum limit being achieved by extrapolating K — > oo. We found it particularly useful 
to employ a Lanczos algorithm in computing and diagonalising the mass matrix M 2 in this way. 
A Hilbert space of dimension ~ O(10 3 ) can be comfortably handled. The elementary processes 
(fig. 1) for adjoint particles represented in (^) include the linear Coulomb potential (E integral), 
present also for fundamental representation particles but with doubled strength g 2 — > 2g 2 since an 
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adjoint source has two flux lines attached to it rather than one. The glueballs (||) may therefore 
be pictured as a ring of flux created by a closed chain of gluons. In addition in ([?]) there is a 2 —* 2 
annihilation channel {D integral) and pair creation and annihilation of gluons (F integral), which 
are not suppressed by the large- N limit but for low- lying levels are kinematically suppressed Q. 
A survey of the broad features of the mass spectrum M 2 resulting from adjoint scalars in 2D has 
been performed in refs.[]l], J3|. A massive spectrum of stable glueballs is organised into approximate 
valence gluon trajectories for low-lying levels, while at higher energies pair creation dominates to 
give a complicated picture. We have concentrated on the highly structured eight or so lowest states 
resulting from ([?]) in order to compare with Monte Carlo results of the unreduced 3-dimensional 
lattice theory. In addition to measuring the masses M from extrapolation of finite K calculations 
to the continuum limit K = oo, we also used the light-front wavefunctions to compute structure 
functions of these mass eigenstates at fixed cut-off. In particular the quantities 

oo „! / n \ 

<n> = Vn/ dxx . . .dx n 5 Vx 4 - 1 |/„(xi, . . . ,x n )\ 2 , (11) 

co „i n 

d( x ) = X] / dxi . . .dx n 5 (j^Xi - l)^25(xi - x)\f n \ 2 (12) 

n=2 J ° V i=l 

give the average number of gluons and the number of gluons with momentum fraction between x 
and x + dx in a boundstate. These quantities are helpful in classifying states according to expected 
quantum numbers of the unreduced 2 + 1-dimensional theory. 

3 Solutions. 

The numerical solutions for the reduced theory are given in table 1 and illustrated on fig. 2 and 
fig. 3. As in refs.Q, [3| we note that a spectacularly accurate valence gluon approximation is at 
work in the light states. The Coulomb potential, which includes the 2nd term in (|7|) as well as 
the ^/-integral, dominates since it is the only positive definite and singular amplitude. Both the 
annihilation channel (D) and pair production (F) amplitudes take either sign with roughly equal 
probability, leading to much destructive interference. Physically they are suppressed because the 
intermediate A+ 'particle' is non-propagating in light-cone time x + , so couples unfavourably to 
physical transverse gluons pairs in the boundstate. 

In labeling the states with their (Pi, C) quantum numbers, a technical problem arises in de- 
termining Pi, which has the effect x + «-> x~ on vectors. This symmetry is broken when the 
momentum fractions x are discretised, since this is equivalent to (anti)periodic conditions on x~ . 
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Under Pi 



1 



(13) 




which in general is not of the form integer /K. There are various quantitative and qualitative ways 
of assessing by inspection the Pi of mass eigenstate at fixed cut-off K. The qualitative one we 
found particularly useful was to note that in general Pi transforms Fock states where one gluon 
carries most of the momentum to Fock states where the momentum is shared evenly between all 
the gluons. 

A striking property of fig. 2 is the almost linear M versus < n > trajectory [|, which is 
approximately repeated before becoming diluted by pair production effects at higher mass. Such 
a trajectory represents 'radial excitation' of the glueball flux loop in the sense that the mass of 
a state is increased by adding gluons and their attendant flux lines to the ring. The leading 
radial trajectory shows featureless structure functions (fig. 3) peaked at 1/ < n >, while the higher 
trajectory exhibits further oscillatory behaviour analogous to higher 'angular momentum' states 
(therefore the Regge trajectories would run vertically in fig. 2). For the latter there is an increased 
probability to find an asymmetrical sharing of the momentum between gluons in the glueball. It 
is interesting to note that g{x) remains large for x ~ in general. In particular the approximately 
two-gluon states are rather like cosine wavef unctions, which was also found for the low- lying mesons 
in the 't Hooft model at small quark mass @. In the latter case the massless groundstate meson 
formed from massless quark and anti-quark had a constant wavefunction f2(x, 1 — x) = const., 
corresponding to non-interacting quarks. In the present case the approximate (++) glueball 
groundstate wavefunction f2(x, 1 — x) = g{x) is not quite constant due to the extra spin-related 
contribution 7r/4y / xiX2 to its mass in (0), suppressing it at the endpoints; this gives the glueball 
spectrum its mass gap. Indeed a pretty good estimate of the groundstate mass is afforded by first- 
order Rayleigh-Schrodinger perturbation theory with f2(x, 1—x) = 1 as zeroth-order wavefunction, 



yielding M ~ 2.2y/Ng 2 /tt, differing from the numerical solution by 5%. In fact all the wave- 
functions on the leading radial trajectory are well described by free gluons f n (xi,X2, ■ ■ ■ , x n ) = 1, 
yielding g(x) = 2n(l — x), plus suppression at the endpoints due to the 1 + 1-dimensional remnant 
of spin interactions, which tend to favour a symmetrical sharing of the momentum. 
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■ 2 Ntt 



(14) 
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4 Comparison with Lattice Results. 

The most reliable lattice Monte Carlo data for SU(3) glueball masses in three dimensions || is 
shown on fig. 4. Less reliable data exists for states nominally of higher mass than those shown. The 
classification is J , where parity P in two space dimensions is taken to be reflection about one 
axis (e.g. x 1 — > —x , x 2 — ► x 2 ). Particles of non-zero angular momentum J should be degenerate 
in P-doublets | J > ±| — J >. Although we have no way of telling the true strength of breaking 
of SO (2) rotational symmetry by the dimensionally reduced theory, which reduces it to a Z2 
subgroup, we will attempt to match our states with J PC labels. Firstly we take P = P\. Under 
180° rotations x 1 — » — x 1 and x 2 — > —x 2 , and although <p is a 1 + 1-dimensional scalar, it is the 
component of a 2 + 1-dimensional vector, so we will assume this induces P\ and <j) — > —(f). Since 
all the states we study are found to be invariant under 4>ij — > <j)ji, which reverses the ordering of 
gluons around a flux ring, this implies that CP\ = \J\ mod2, which rules out half the possible 
states in three dimensions from the very beginning. In fact it removes precisely one of the states 
from each P-doublet. One is still left with the problem of spin labeling within the set of even and 
odd spins separately, because in general there will be mixing within each set, i.e. spin mixes 
with spin 2 etc.. We took the spin label for states of given (Pi, C) on a radial trajectory to be the 
lowest which has not appeared on a lower energy trajectory with the same (Pi, C). 

The mass ratio scale gets expanded in the truncated theory, presumably due to the mixing 
between levels as a result of 50(2)-breaking, so we have shrunk it by an appropriate factor in fig. 4 
in order to better compare the level ordering with lattice results. It is tempting, in this gauge 
and in this large- iV light-front quantisation scheme, to describe the data in terms of a constituent 
gluon picture. However, in at least one case this is in clear conflict with the {N = 3) lattice 
data: a O state could be interpreted as an orbital excitation of the 3-gluon groundstate l" 1 ; 
on the lattice the O is the lower state since it couples to the combination Up — U P of one 
Wilson plaquette Up, while l" 1 needs a longer lattice Wilson loop. In fact fig. 4 shows that the 
(N = 00) dimensionally reduced spectra in the C = — 1 sector predict light spin 1 glueballs, in 
disagreement with the 3D lattice. The ordering in the C = +1 sector is quite good by comparison. 
Spin l's being lower than spin 2's may be due to the mixing of even spins and odd spins amongst 
themselves since spin 1 is pushed down by mixing with spin 3 etc. 

Lattice results also exist for SU{2) when C = +1 []|], but states with C = — 1 do not appear 
since Wilson loops are unoriented. For completeness we have therefore plotted an N = 00 extrap- 
olation for C = +1 from N = 2,3 data by fitting to the form = Mn — const. /N 2 , though 
the variation is actually quite small. We do not think that finite N corrections are a major source 
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of error, relatively speaking, but that the discrepancies we have found are probably largely of a 
kinematic nature due to the severity of the reduction and hence the smallness of the residual Z2 
symmetry. We are however encouraged by the simple physical picture that the reduced theory 
presents since it incorporates the two principal dynamical effects at work in hadrons, the linear 
string potential and spin-dependant interaction. Although the reduced model also includes pair 
production of transverse gluons, the explicit computation shows this to be suppressed to a re- 
markably high degree. There is no direct self-coupling of transverse gluons in three dimensions, 
only coupling through a longitudinally polarized non-propagating gluon, which we believe is the 
origin of the suppression of pair production. The situation when reducing gauge theory from four 
to two dimensions is rather different and it is clearly of interest to perform a similar analysis for 
the 1 + 1-dimensional two-adjoint scalar theory obtained by ignoring XT-dependence. It has the 
light-front hamiltonian (x^ = (x° ± x 3 )/\/2) 



where fa ~ Ai are the transverse zero modes of the gauge potential, and includes a 4-point contact 
interaction between transverse fields. The discrete remnant of the rotation group in this case is 
the much larger Dihedral group D4, so there is some hope that the mixing between different spins 
is lessened. A detailed analysis of the hadron spectrum, including quarks, will be given elsewhere. 
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FIGURE CAPTIONS 
Figure 1 - The four elementary processes contributing at order g 2 in Sr. 

Figure 2 - Mass spectrum extrapolated to the continuum limit K = oo, classified by (Pi, C). The 
estimated error in these masses is a few percent, though great accurracy is not terribly important 
for the purposes of this paper. Note that the P\ eigenvalue of the highest state shown could not 
be identified with certainty (see text). 

Figure 3 - Structure functions of mass eigenstates: (a) < n >~ 2 ; (b) < n >~ 3 ; (c) < n >~ 4 ; 
(d) < n >~ 5. The solid lines refer to states on the lower radial trajectory in fig.2 while the chain 
lines refer to states on the upper trajectory. 

Figure 4 - Comparison of dimensionally reduced and 3D lattice glueball mass ratios to the ground- 
state. SU (3) lattice data has solid circles, the N — > oo extrapolation of lattice data where available 
has open circles, and the reduced theory has squares. Note the different scales; that in parentheses 
refers to the reduced theory results. 

Tablel - M is the K = oo extrapolated mass using the full basis of states. < n > is calculated 
from the theory truncated to the sector of 2, 4, and 6 gluons only for K = 20, or to the sector of 
3 and 5 gluons for K = 19. 
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